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A putative powerlaw range of the probability density of velocity gradient in high-Reynolds-numbcr 
forced Burgers turbulence is studied. In the absence of information about shock locations, ele- 
mentary conservation and stationarity relations imply that the exponent —a in this range satisfies 
a > 3, if dissipation within the powerlaw range is due to isolated shocks. A generalized model 
of shock birth and growth implies a — 7/2 if initial data and forcing are spatially homogeneous 
and obey Gaussian statistics. Arbitrary values a > 3 can be realized by suitably constructed 
homogeneous, non-Gaussian initial data and forcing. 



Burgers equation was originally proposed as a sim- 
plified dynamical system that might point at statistical 
procedures applicable to Navier-Stokes turbulence. That 
goal seems still far away. Meanwhile, a body of research 
has been devoted to Burgers turbulence itself. The par- 
ticular topic addressed in the present paper is the power 
law of a part of a probability density for Burgers turbu- 
lence forced at high Reynolds number. The conclusion 
reached is that determination of the correct power law 
requires detailed examination of the dynamical behavior 
of explicit structures, namely shocks. This perhaps has 
discouraging implications for any general theory of the 
higher statistics of Navier-Stokes turbulence, where more 
varied and plastic flow structures must be confronted. 

There has been considerable interest in a putative pow- 
erlaw range of the probability density Q(£) of veloc- 
ity gradient £ = u x for high-Reynolds-number Burgers 
turbulence forced at large scales. This range, of form 
Q(£) 00 l£l -Q S i s believed to occupy negative £ values in- 
termediate between those near the central peak of Q and 
those characteristic of the shock interiors. Proposals for 
a include 2 @J|, the range 5/2 to 3 3 §, and 7/2 

§§■ 

It should be emphasized at the outset that the high- 
Reynolds-number limit is not the only case of theoreti- 
cal interest. As with Navier-Stokes turbulence, the con- 
struction of faithful analytical approximations at finite 
Reynolds numbers remains a challenge. 

E and Vanden Eijnden advance and clarify the 

mathematics of the infinite-Reynolds-number limit and 
provide some valuable tools. One is a simple represen- 
tation of the effects of shock interactions on Q in terms 
of the rate at which fluid is swallowed by the shocks. 
Another is a steady-state integral representation of the 
asymptotic large- 1£| form of Q(£) in terms of the dissipa- 
tion term F(£) in the Q equation of motion. Another is 
an explicit shock-birth model that implies a = 7/2. 

The present paper explores the relation between the 
form of F(£) and the value of a. If viscous effects in the 
—a range are due to isolated shocks, the form of F(£) in 



this range expresses the relative likelihood, weighted by 
shock strength, that a shock occurs in a fluid environ- 
ment with given £. 

In the absence of an explicit shock-growth model, or 
other source of information about the distribution of 
shocks, the integral equation for Q is found to yield 
a > 3. The slightly stronger bound a > 3 is stated 
in and g, but with the recognition that a — 3 may 
be realized under particular circumstances. The limit of 
infinite Reynolds number is taken in the present paper 
without making a split jl]] of the velocity field into shock 
interiors and external field. At the end, the analysis is 
extended to the split-field representation. 

A more general model of shock growth is presented 
here. It is independent of details of internal shock struc- 
ture. In this model, one examines the length of time 
during which £ can steepen within a fluid element be- 
fore the fluid element hits a shock. The model implies 
a = 7/2 if the forcing is statistically homogeneous and 
Gaussian. More general statistically homogeneous forc- 
ing can realize arbitrary values a > 3. 

Let R — uqL/v, £o = uq/L, £5 = i?£o, where uo, L, 
v are root-mean-square velocity, spatial macroscale, and 
viscosity. The order of magnitude of gradients within 
typical shocks is £g. The forced Burgers equation is 

u t + u x u = vu xx + /, (1) 

where / is the (statistically stationary) forcing field. If / 
has infinitely short correlation times, (]l|) leads to 

Qt = £Q + (£ 2 Q)c + BQ 55 + F, (2) 

where 

F(t,t) = -v(H{t,t)Q(t,t)h, = <6«I0> (3) 

and the parameter B measures the strength of forcing of 
£ [pi. The first term on the right side of (||) represents 
loss or gain of measure due to squeezing or stretching of 
the fluid, the second term describes relaxation of positive 
£ and steepening of negative £, and F includes all viscous 
effects. 
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Statistical homogeneity requires (£) = J™ £,Qd£ = 0. 
Multiplication of (||) by £ and integration over all £ shows 
that this condition is preserved, provided that Q vanishes 
strongly enough at ±00. The result depends on 

/oo 
£F(£)d£ = 0, (4) 
-oc 

which follows from (||) and (£ X x) = 0. 

If / has spectral support effectively confined to 
wavenumbers 0(1/L), then £0 = 0(B 1 ^ 3 ). In this case, 
it is widely agreed that the steady-state Q has a com- 
plicated form in the limit R — > 00. There is a central 
peak of width O(£o) whose form is R- independent in the 
limit. There is faster-than- algebraic decay as |£| — > 00. 
For £ > 0, this decay has the specific form exp(— £ 3 /3B), 
with an algebraic prefactor. The rapidly-decaying tail 
for £ < (far- left tail) includes |£| > 0(£s). It is pre- 
ceded, at smaller |£|, by an algebraic tail of form 
(— 1 range) associated with the shoulders of developed 
shocks. Between the —1 range and the central peak, an 
inner algebraic tail of form Q(£) oc £g _1 |£|~ Q is expected. 
This tail is driven by the inviscid steepening of negative 
gradients. Proposals for the value of a have ranged from 
2 to 7/2. 

The —a range is infinite if R is, but is confined to 
|£| smaller than the 0(£s) gradients inside the shocks. 
Thus the range is restricted to £ > — £m = ~R z £,o, with 
z < 1. In fact, both shock analysis and simulation show 
that the —a range is masked by the —1 range at negative 
enough £. The masking further restricts the observable 
—a range to z — l/(a — 1), a result that follows from 
setting 1/R£m = Co" 1 ^"- ^ snou ki be emphasized 
that the transition from Q(£) oc |£|~ Q to oc 
at £ = 0(— £m) is a masking, not a dynamical transi- 
tion. The —I range is associated with the shoulders of 
quasi-stationary mature shocks while the —a range is as- 
sociated with inviscid steepening of gradients away from 
shocks. The latter process continues, for at least some 
fluid elements, up to |£| = 0(£s). 

The behavior of Q in all ranges is linked by (||) to the 
form of F(£). For large negative £, the relationship is 
expressed in especially clear form by the integral repre- 
sentation 

Q(o*fcr 3 / mew (-£»&) w 
j — 00 

derived by E and Vanden Eijnden |Q. An alternative 
form of (§) is || 

^(0«£ 2 + ^y (-£»&)• (6) 
or by (|), 

^«-£Q-(£ 2 Q)e (-£»&)• (7) 



Equation (|7j) says that in a steady state, for £ such that 
the B term is negligible, the F term in (||) must balance 
both the gradient steepening term (£, 2 Q)^ and the loss- 
of-measure term £Q. Equations (^)-(Q) do not require 
large R. 

According to (|1), the behavior of Q in the —a range 
depends on, first, the value /_ M £'F(£')d£' of the inte- 
gral at the outer end of the range; second, the £ depen- 
dence of F(£) within the —a range, which determines the 
growth of the integral within the range. If F(£) in the 
—a range arises solely from isolated shocks, the analysis 
of |Q] gives an immediate bound on a. This analysis ex- 
presses a simple physics: If a shock lies in surrounding 
fluid with gradient £, then this fluid is swallowed by the 
convergence at the shock, thereby decreasing Q(£). This 
means that F(£) cannot be positive in the —a range. 
Since £ < 0, the integral in (@) therefore cannot increase 
as |£| decreases within the range, and a < 3 is impossible. 

The cases a = 3 and a > 3 imply qualitatively differ- 
ent magnitudes in (||). If a = 3, then (0) immediately 
gives XT^ M £'F(£')d£' = O(fg). It is then required that 

F be small enough that £'F(£')d? = 0(g) if £ lies 
anywhere within the —a range. If a > 3, (g|) immedi- 
ately gives J^' t'Fffld? - as £ M /£o -» 00. Then 
F(£) oc — |£| 1_Q within the range will give Q(£) oc 
This means F(£) oc £Q(£)- 

If dissipation within the —a range is due to isolated 
shocks, the physical difference between a = 3 and a > 3 
lies in the weighted relative likelihood that shocks exist 
in environments with given £. -F(£) oc — <2(£), a case 
of Polyakov's closure ||,|} , is one form that is consistent 
with a = 3. This corresponds to unbiased placement of 
shocks in the fluid jnj. -F(£) oc £,Q(0 represents bias 
toward larger |£| within the —a range. Since Q(£) falls 
off rapidly in the —a range, even this form puts most of 
the shocks in fluid with |£| = O(£o)- 

If a = 3, the inviscid steepening of gradient in most 
fluid elements with negative £ continues until |£| = 
0(£s). If the dissipative process that stops further steep- 
ening is absorption by well-defined shocks, these shocks 
must occur with sufficient weight in environments with 
|£| = 0(£s)- Since £s measures typical gradients within 
shocks, other processes, such as new shock formation, 
may actually dominate the dissipation. 

Further insight concerning the relation between F and 
a is provided by the equation of motion (||) for Q and, 
in particular, by the steady-state equation Q t — 0. A 
variety of global assumed forms for F, including F(£) oc 
— Q(£), can be adjusted, by tuning of the constant of 
proportionality, to give a steady-state Q(£) that vanishes 
strongly as £ — > +00, is positive everywhere, and has a 
range a = 3. In order to be physically relevant, F(£) for 
negative £ beyond £m must be shaped to give consistent 
1/-R|£| and far-tail ranges for Q. Global model forms can 
also be constructed that yield a range with a > 3. 
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The quantity (£) is identically conserved at value zero. 
It is of interest to examine the way in which flow of 
(£) through the —a range depends on the value of a. 
Let (^|) be multiplied by £ and integrated over the range 
(-£ M ,oo) to yield 



£Q t (£R = e M Q{- 



£m) + 



SF(S)d£. (8) 



In writing (g), a partial integration is performed, it is 
assumed that Q vanishes strongly enough at +00, and it 
is assumed that £m/£o is large enough that the B term 
is negligible. 

The left side of (||) is the rate of increase of (£} in 
the range (— £m>°o)- On the right side, ZmQ(~Zm) is 
the rate of increase due to flow of negative contribu- 
tion to (£), through the boundary at — to the range 
(—00, — £m)- This flow is due to inviscid steepening of 
gradients. The F term on the right side of (||) is the 
rate of increase of (£), or decrease of (— £), in the range 
(-(if, 00), due to viscous interaction with isolated shocks 
and any other dissipative structures that may be present. 

The nature of the F term in (||) must be understood 
clearly. The present analysis is at finite R, with the even- 
tual limit R — > 00 considered. There is no sudden jump of 
(£) contribution from —a range to shock interiors. What 
does happen at large R is that a fluid element with £ 
in the —a range hits the shock and then suffers a very 
rapid, but continuous, steepening until its £ is the order 
of that in the shock interior. The sum of these events in 
the entire range is expressed by the F term in (^). 

In a steady state, the right side of (||) vanishes. The 
implications differ in the cases a = 3 and a > 3. If 
a > 3, the boundary- flow term in (^) vanishes in the 
limit £m/£o — * °o- This means that the F integral term 
vanishes also in the limit. In view of where the lim- 
its are true infinity (infinite compared to £g), it follows 
that /_ £F(£)dS; also vanishes in the limit. This does 
not mean that F(£) tends to zero in the limit for all £ in 
the range (— 00, — £m)- Instead, there are both positive 
and negative contributions that cancel in the limit. In 
general, F (£) is positive in the —1 range, as illustrated 
by © an d djjj ) below. 

If a = 3, the boundary term in (^) tends to a nonzero 
positive constant in the limit. A steady state then im- 
plies that the F integral in (||) is negative and that 

f-t* Z F (£) d Z is positive. This means that the contri- 
butions from the negative and positive regions of F in 
the range (—00, — £m) do not cancel. The nonzero value 
of j:^ 1 £F(£)d£ needed in the case a = 3 has already 
been noted in the discussion of (|J). 

For both a > 3 and a = 3, the boundary flow through 
an arbitrary point — £ Q fully within the —a range is in- 
dependent of £ Q in the limit R — > 00. Thus if £a/ in (^) 
is replaced by any £ Q such that £ q /£a/ — ► 0, £ Q /£o — ► 00 
in the limit, then the limiting value of the boundary flow 



vanishes if a > 3. The flow has a value 0(£ 2 ), indepen- 
dent of £ Q , if a = 3. 

The form of Q and F in the far- left tail and —1 range 
actually is not constrained by whether J_ £-F(£)d£ 
vanishes or is O(£o) as R — > 00. The key is the value 
of z = l/(at — 1), which gives the position of the join 
between —1 and —a ranges at — £m = ~R z Zo- Consider 
the generic form 



Q(£) w Z(£/i?£ )/i?|£| (-£ » £ ) 



(9) 



in the far-tail and —1 ranges, where i? vanishes strongly 
as £ — > — 00, Z — > C 2 as £/i?£o — > —0, and C 2 is an O(l) 
constant. If this form is substituted into (0), a conse- 
quence is 
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£F(£)d£ « C 2 i?^£ 2 (i?»l) 



(10) 



where /3 = (3 — Q>)/(a — 1). If a = 3, this gives the 
needed O(£o) result. If a > 3, the right side vanishes at 
R = 00. Thus the value of the boundary flow automati- 
cally adjusts to the value of a that is determined by the 
form of F in the —a range. In other words, it adjusts 
to the probability distribution of shock occurence in the 
—a range. 

Explicit models of shock development lead to explicit 
forms of .F(£). The following model generalizes those of 
. It assumes large R but does not invoke the internal 
shock structure. The development is followed before and 
after shock birth. Take the unforced case first. Consider 
an initial velocity field of form 



u(x,0) K^ (~ax + bx\x/L\ p ) (p > 0) 



(11) 



in the vicinity of a point x — where a shock will form at 
time l/a£o- Here a and b are O(l) positive constants. If 
p is an even integer, all x derivatives of u exist at x — 0. 
If p is an odd integer or non-integer, only the derivatives 
of order n < p + 1 exist. 

This form of initial velocity field leads to a — 3 + 1 /p, 
a result that can be verified in several ways. The fol- 
lowing is a simple qualitative argument. The inviscid 
evolution of velocity gradient in a Lagrangian frame sat- 
isfies £ = — £ 2 . It then follows from ([ll]) that the time 
of initial shock formation is to = l/a£o- The negative 
gradient within a fluid element initially at small enough 
\x\ > grows inviscidly until the fluid element hits the 
shock. 

The time at which a fluid element initially at x falls 
into the growing shock is t x ~ l/[£o(a — b\x/L\ p )}. The 
times of arrival at the shock determine the fractional 
measure of initial points x such that, before a fluid el- 
ement hits the shock, the gradient magnitude increases 
to values > |£| 3> £o- The gradient magnitude at x grows 
to w a£ \x/L\-P/[b(p+ 1)] at t and w a£ \x/L\-P/bp at 
t x . Thus the measure of points x such that the gradient 
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magnitude in the fluid element initially at x can grow to a 
value that equals or exceeds |£| 3> £o within the intervals 
(0, to) or (0, t x ) is oc l^l" 1 ^. Of the initial fluid elements 
that achieve a value at least |£| 3> £o before hitting the 
shock, the fraction that does this in the preshock interval 
(0,t ) is ( P /p+l) 1/p . 

When the fluid element has achieved the value £, 
squeezing has decreased its measure by a factor oc 
Finally, the residence time of the fluid element in the 
gradient interval d(, at £ is dt = |£| _2 d£. Putting these 
factors together, one obtains Q(£) oc \£\- 1 - 2 ~ 1 /p = 
|£|~ 3 ~ 1 / p , where £}(£) is the mean of Q(£, t) over a time 
interval (say 2/a£o) long enough for all fluid elements 
that can achieve |£| ^> £o to have hit the shock. 

In ||, ol = 3 was deduced under the assumption 
that the fractional measure of fluid elements that can 
achieve gradient magnitudes » £o is 0(1). The measure 
oc |£| _1 / p found instead in the present model changes the 
result to a = 3 + 1/p. For all finite p > 0, the form of 
F calculated from the present model is F(£) oc 
within the —a range. 

Initial sawtooth profiles, where u(x) consists solely of 
straight-line segments, correspond to b = (or p = oo) 
in the model. They evolve into shocks that have finite 
amplitude at birth and yield a — 3. The present model 
thereby is consistent with the conclusion that iso- 
lated shocks can induce a = 3 only if they are created 
with finite amplitude. 

The steady state produced by spatially smooth Gaus- 
sian forcing supported by wavenumbers 0(1/ L) can be 
interpreted in terms of this model. Such forcing induces 
smooth profiles corresponding to p — 2 near points of 
extremal slope. In the absence of force, a shock forms 
from steepening of slope at a point of maximally nega- 
tive slope. Smooth change of the velocity field due to 
Gaussian forcing can change the location of the point of 
maximally negative slope as a function of time. But such 
forcing does not change the nature of the shock forma- 
tion phenomenon because the quadratic decrease of slope 
magnitude away from the point of negative maximum 
survives. Once the slope magnitude at negative maxi- 
mum is large compared to t;o, the forcing should have no 
significant effect on either the progression to shock birth 
or the initial shock growth after birth. The value a = 7/2 
corresponding to p = 2 is the steady-state result. 

Forcing that supports the general case p ^ 2 in a sta- 
tistically steady state can be constructed as follows: Let 
the forcing consist of a set of ^-functions in time, spaced 
at time intervals O(l/£o). Let each such 5-function cre- 
ate an increment to u(x) that consists of straight-line 
segments of length 0(L) with O(u /L) positive slope, 
smoothly joined to interposed negative-slope regions of 
O(L) or shorter lengths. In each negative-slope region 
let there be a point of maximally negative slope sur- 
rounded by a neighborhood in which u(x) has the form 
£,o(— ay+by\y/L\ p ), where y is the distance from the point 



of maximally negative slope. The values of a and b can 
change stochastically from one such region to another. 
Under these conditions, the negative-slope increment to 
u(x) created by each 6- function force field is added to an 
existing field that has locally constant slope with 0(1) 
probability. Therefore the points of maximally-negative 
slope are at the special points y — 0, and the consequent 
shock development yields a = 3 + 1/p for all p > 0, as in 
the initial-value case. 

All cases of the model except p — 2, a — 7/2 require 
special shapes of the velocity field prior to shock for- 
mation and, therefore, precise phase relations of spatial 
Fourier components. If the forcing field is spatially ho- 
mogeneous and has an infinitely short coherence time 
(white forcing), the effective forcing is Gaussian and 
these shapes and phase relations cannot be realized. The 
white, homogeneous forcing assumed in writing the B 
term in (0) implies p — 2, a — 7/2. 

At the end of Q , it was noted that the value of a de- 
pends on how likely it is for a shock collision to interrupt 
the steepening of negative gradients. It was argued that 
collisions should be infrequent enough that inviscid steep- 
ening should survive for most fluid elements of negative 
£ until \£\ = 0(£s) is reached. This implies a = 3. The 
present shock-growth model, following the earlier ones in 
HQ, says instead that fluid elements with large nega- 
tive £ are inevitably close to shocks that soon swallow 
them if p = 0(1). For most fluid elements with negative 
£, the inviscid steepening is terminated by collision with 
the shocks while |£| is still O(£o). 

The analysis in |Q is done in terms of a split of u and £ 
into parts exterior to shocks and parts interior to shocks, 
in the limit R — > oo: 

u(x, t) = u e (x, t) 4- Ui(x, t), £(x, t) = £ e (x, t) + t) 

(12) 

It is of interest to discuss how to make the split (|l2| ) when 
R is large but not infinite, and to express the preceding 
analysis in terms of the split-field representation. The 
split into interior and exterior fields is analyzed in |||]] 
by means of matched asymptotic expansions. 

If the R — oo field consists of infinitely thin shocks 
surrounded by fluid in which |£| < £i, where £i is some 
finite bound, the split is clear and unambiguous. At large 
finite R, one can isolate a small region surrounding each 
shock [Q as Ui, and consistently let the widths of these 
regions shrink to zero at R — oo. E and Vanden Eijnden 
show that the result is a simple set of statistical equations 
relating £ e and the shock jumps: 

(&>+P (*>=(>, (13) 
F*&,t) = % I s[V^ ei s,t) + V+(^ S ,t)}d S . (14) 
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Here p is the number density of shocks, s is shock-jump 
strength (< 0), F e is the viscous term in an equation of 
motion like (||) for the probability density Q e (£ e ) of the 
exterior field, and V- (V+) is the probability that a shock 
of strength s has a left (right) environment with gradient 

These equations have a direct physical interpretation. 
Equation ( |l4] ) expresses the viscous term as the rate at 
which Q e (£e) is diminished through the swallowing of 
fluid with gradient £ e by shocks. The time derivative of 
( |l3| ) is an expression of the fact that the rate of change of 
shock strength is given by the product of the convergence 
velocity and the negative of the gradient of the external 
field, as the latter is swallowed by the shock. 

The split into £ e and £j is less clear if there is a —a 
range at R = oo that includes |£| values larger than 
any finite bound. Then the £ e field is not smooth. As 
R approaches infinity, it is not possible to form an in- 
dependent parameter £1 such that all field with |£| < £1 
belongs to £ e and all field with |£| > £1 belongs to £j. 
Instead the division into interior and exterior fields must 
be made individually at each shock. 

The equation of motion for £ obtained by differentia- 
tion of (Q) contains the steepening term £ 2 and the dis- 
sipation term v£ xx . One possibility for splitting the £ 
field is to let £ e include all points £ > and all points 
£ < that are exterior to defined boundaries of the shock 
shoulders. A shock-shoulder boundary could be defined 
as a point where, as one moves through the field toward 
the shock, t; 2 /\v£, X x\ first falls below some prescribed ra- 
tio, say 10. Then & constitutes the field interior to the 
left and right shock-shoulder boundaries. This gives a 
£ e field that extends up to values |£| = 0(£s) at some 
points, with a Q e (£e) that must deviate from power law 
behavior for such |£| values. The qualitative behaviors of 
£ e and £j are independent of the precise value prescribed 
for the critical ratio of steepening term to dissipation 
term. 

Since |£ e | extends to oo in the limit, it is not obvi- 
ous that all interactions between exterior and interior 
fields have a form consistent with (jlj) and In the 

limit, the 1/-R|£| range belongs to which includes the 
shock shoulders. Both £ e and contribute to the total- 
field probability density Q(£m), and this fact remains as 
R — > oo. 

It is been remarked above that F(£) is positive in the 
— 1 range, while (14) gives negative F e (£ e ). This empha- 
sises that the region of £ contributing to the — 1 range in 
each individual shock boundary layer should be assigned 
to 

Equations similar to (g), (||), and (||) can be written 
for Q e (£ e ): 



CF e (Q< K»(o), (16) 



(17) 



Although the equations look the same, F e and F behave 
differently for negative arguments. 

The parameter £m has no special significance here be- 
cause the —a range of Q e is not masked by the —1 range. 
The latter belongs to . Therefore £m has been replaced 
in (^) by £ a , which is any value within the —a range 
that satisfies the limiting relations 



° (iZ->oo) 



(18) 



If a = 3, (|17|), like (g), exhibits a boundary flow that 
does not vanish at R = oo. 

The analog of (g), £ e F e (( e )d£ e = 0, holds in 
steady states. In general transient states, (|l3| ) implies 
that there is an additional term involving the rate of 
change of shock jumps. 

The power-law behavior of Q e {S,e) must change to a 
faster (eventually faster-than-algebraic) fall-off for |£ e | > 
0(£s). The generic form may be written 

QeiQ-Hg-^el^Z^/RZo) (-&»$>)• (19) 

Here Z vanishes strongly at £ e = — oo and is 0(1) at 
£ e = —0. The precise form of Z is oj-dependent. In con- 
trast to (^), there is no —1 range. The prefactor in ( |l9| ) 
gives the consistency property Q e {— Co) = O(l/£o) if tric 
—a range is extrapolated toward the central peak of Q e . 
Substitution of dH) into © yields 



"Co 



(20) 



(Qe)t = LQe + (&Qeh + ^(Q e )« + F e , 



(15) 



If a > 3, the dissipation measured by £ e -F e (£ e )(i£ e 
equals the total dissipation of £ e in the limit R — > oo. 
If a = 3, ( p0| ) shows that there is additionally an essen- 
tial contribution J_ £ e F e (£ e )d£ e that is 0(£q) in the 
limit. This arises from 0(£o/£s) levels of -F e (£e) needed 
to induce the fast fall-off of Q e (£e) at |f e | = 0(&?)- The 
fast fall off occurs also if a > 3, but in that case J2(i| ) 

shows that the associated contribution S, e F e {^e)dS, e 
is a vanishing part of the overall budget in the R — > oo 
limit. 

In [Q, the limit "— oo" in ( |l6| ) above and other equa- 
tions is taken to mean a point within the infinite —a 
range, rather than true negative infinity (<C — £s). In 
other words "—oo" is taken to mean — £ a , as defined by 
( |l8| ) in the limit R — oo. This causes no problem if a > 3. 
If a = 3, it is clear from ( pp[ ) that "— oo" must stay at 
true negative infinity. 
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If the possibility a = 3 is to be examined, clearly one 
cannot confine attention solely to the strict —a powerlaw 
range and ignore the transition region at |£ e | = 0(£s). 
The case a — 3 implies that, for most fluid elements with 
steepening negative £, the inviscid steepening halts at f; 
within the transition region. The magnitude of -F^^e) at 
|£ e | = 0(£s) required by a = 3 signifies that shocks have 
a relatively high concentration in fluid with such £ e . In 
comparison with a > 3, shocks are moved from environ- 
ments with |£ e | <C £s to environments in the transition 
region |f e | = 0(£s). 

If an explicit shock-growth model is not adopted to fix 
F e , it is possible a priori that shocks with environments 
in the transition region could fall within the picture in- 
voked by |t]], in which shocks are created at zero ampli- 
tude and the balance is described fully by ( |l^ ) and (14). 
In this case it is only needed to include the transition 
region of shock environments in (£ e ) when calculating 
the balance ( |l3| ) for a = 3. However, it is also possible 
a priori that F e (£ e ) is weighted to sufficiently large |£ e | 
that ( |l4]) is not an accurate description of the dissipation 
mechanism for a = 3, or that -F e (£e) includes interactions 
with shocks created at finite amplitude as also considered 
in §. 
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